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Abstract 

The study of B —* r + r~ (X) decays can provide us with a better under- 
standing of the third generation, and can be a useful probe of physics beyond 
the standard model. We present a model independent analysis of these de- 
cays. We classify new physics that can largely enhance the decay rates and 
we discuss the constraints implied by other processes. Experimentally, flavor 
changing neutral current B decays into final state r's are still unconstrained. 
Searches for B decays with large missing energy at LEP provide the first lim- 
its. We estimate that existing data already imply bounds on the B^ — * t + t~ , 
B s — * t + t~ , and B — > X r + r~ decay rates at the few percent level. Al- 
though these bounds are over four orders of magnitude above the standard 
model predictions, they provide the first constraints on some leptoquarks, and 
on some R-parity violating couplings. 
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I. INTRODUCTION 



The standard model (SM) of the strong and electroweak interactions provides an accu- 
rate description of the low energy properties of the known elementary particles. However, 
experimental results involving fermions of the third generation are far less precise than for 
the first two generations. For example, little is known about decays involving more then 
one third generation fermions. From the theoretical point of view, a better knowledge of 
the physics of the third generation could help us understand the hierarchy of the fermion 
masses and mixing angles. As it is the case in some models of new physics [[!]], the third 
generation might even be essentially different from the first two. 

The LEP experiments have provided us with several new results on third generation 
fermions, and in particular on the b quark. One type of measurement for which the LEP 
environment has advantages over symmetric B factories (like CLEO) or hadron colliders 
(like CDF) is the study of B decay modes that produce a large amount of missing energy 
due to neutrinos in the final state. The main background to such analyses is the tail of the 
semileptonic decay distribution, so decay modes yielding a harder missing energy spectrum 
can be effectively measured or constrained. The excess of events with large missing energy 
measured at LEP was interpreted as the signature of the decay B — > X c tv t followed by 
t^vX M, yielding 

BR(5 ^X c tu) =2.68 ±0.34%. (1.1) 

This is in good agreement with the SM prediction BR(5 — > X c r v) = 2.30 ± 0.25% ||, and 
constrains certain new physics contributions ||. The ALEPH Collaboration also searched 
for events with very large missing energy, E m i SS > 35 GeV 0. The absence of excess events 
over the background yielded the 90% confidence level upper limit on the exclusive leptonic 
decay B — > r v [0 

BR(5 -> tu) < 1.8 x 10~ 3 . (1.2) 

In a recent paper |7j we showed that the same data also imply a bound on the flavor 
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changing B — > X v v decay rate, and we discussed the resulting constraints on several possi- 
ble sources of new physics. Based of the full LEP-I data sample, the ALEPH Collaboration 
has recently announced a preliminary 90% confidence level limit on this decay mode || 

BR(5 -> Xvv) < 7.7 x 10 -4 . (1.3) 

This limit is only one order of magnitude above the SM prediction, and provides important 
constraints on several new physics scenarios [fjj . 

Besides these decay modes, the exclusive B — > r + t~ and inclusive B — > X r + r~ decays 
are also associated with sizable missing energy due to the neutrinos from the r decays. Since 
these processes are presently unconstrained, it is interesting to see whether any useful limit 
can be established by analyzing the LEP data on large missing energy events. In the SM, 
B decays into a pair of charged leptons are highly suppressed. However, certain kinds of 
new physics can enhance these decay rates up to several orders of magnitude above the SM 
predictions. 

In section II we study the B^ s — ► r + t~ decays from a theoretical point of view. Since we 
are mainly interested in possible new physics contributions, we present a model independent 
analysis. In section III we estimate the limits on the B — > t + t~ (X) branching ratios that 
could be established from the LEP data. In section IV we discuss the constraints that the 
limits on these decays imply on some new physics models. Section V contains a summary 
and our conclusions. 

II. THEORETICAL ANALYSIS 

The most general effective four-fermion interaction involving a b quark, a q = d or s 
quark, and a pair of t + t~ leptons can be written in the form 

c* =g f y: ci (qr a b) (f r a t) + g f y: cf (q r a b) (t r' a r ) , (2.1) 

a a 

where T a = {1,75,7^,7^75,^}, = r a75 and a = {S,P,V,A,T}. In we have 
factored out the Fermi constant so that the coefficients C% and C%' are dimensionless. 
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While the B — > X q r + t~ decay depends on all ten coefficients C% and C%', fewer operators 
contribute to the leptonic B q — > r + r~ decay.f] 

Let us consider the matrix element (0\qT a b\B). It can only depend on the four- 
momentum of the B meson, p B , and therefore it must vanish for Ft = which is antisym- 
metric in the Lorentz indices. The matrix elements of the parity-even operators (Lg = / and 
Ty = 7 M ) also vanish due to the pseudoscalar nature of the B meson. Finally, for on-shell 
leptons, the contribution of the axial- vector operator (0| #7^75 b \ B) oc p B = p^ + + p^- also 
vanishes when contracted with the leptonic vector current f 7 m t. Hence, leptonic B decays 
are induced only by the following three operators 



C q P (?75 6) (r 75 t) , C q P ' (q 75 b) (f r) , G\ (q 7^75 6) (r W r) . 



(2.2) 



Using the PCAC relations 



(0\q 1 ^ 5 b\B) = -zf B p» B , 



q^b\B) = if B 



m B 
m\ + nig 



ifB rriB , 



the most general amplitude for the B q — > r + r decay reads 



A q = %j B m B G F 



m B 



(2.3) 



(2.4) 



where, for simplicity, we omitted the subscripts q = d,s for the B mass and decay constant. 



All three operators in ( |2.2| ) appear in the SM. Thus, the general result for the total decay 
rate can be read off from Ref. |jj 

G 2 F f B m% I 



T(B q ^T + r- 



4 m 2 



8tt 



ni 



B 



2 m r ,.„ 2 A 4m 2 



4<\ c ,/ 



(2.5) 



In the SM, CJ, and C P gets contributions from penguin diagrams with physical and 
unphysical neutral scalar exchange, and are suppressed as ~ (mb/mw) 2 - Then, the dominant 
contribution to the decay rate comes from 



*We thank David London for bringing this point to our attention. 
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v 1 V8 7rsin 2 ^ w \ m w J 

where, at leading order, Y(x) = (x/8) [(4 — x)/(l — x) + 3x lnx/ (1 — x) 2 ] fTOH - Including the 
small next-to-leading order correction, the SM result for the branching ratio is Iff!"] 



BR SM (£ S -> r+ r~) = 8.9 x 10" 



230 MeV 



m t [rn t ) 
170 GeV 



3.12 / It r I \ 2 



Vu\\ I T B S 



0.040/ V1.6ps. 



• (2-7) 



Compared to this result, the Bd — > r + r decay rate has an additional suppression of about 
\V td /V ts \ 2 ~ 10- 2 . 

In general, new physics that can induce large contributions to the coefficients in (|2.5|) is 
also likely to enhance the rates for other rare processes. The existing experimental limits 
already imply that in several models, the B q — > r + t~ decays can only occur at rates far 
below the sensitivity achievable at LEP. Therefore, it is useful to classify the contributions 
which are already tightly constrained by other measurements, and the ones which are still 
unconstrained (or only weakly constrained). In doing so, we will avoid referring to any 
specific model, but we will use SU{2) gauge invariance to relate operators contributing to 
the B q — > r + r~ decays to operators which induce other transitions. In the presence of new 
physics, operators which are not manifestly SU{2) invariant can also appear []12|, suppressed 
by inverse powers of some large mass scale. They will not be considered here, since it is 
unlikely that through such contributions the B q — > r + t~ decay rates could get the few order 
of magnitude enhancement required in order to be observable at LEP. 

SU(2) invariants can be built out of four SM fermions by combining four singlets, four 
doublets, or two singlets and two doublets. While for l 4 and l 2 x 2 2 there is only one SU(2) 
invariant, two different SU{2) invariants arise from 2 4 . Taking into account all inequivalent 
permutations of the b, q, and r fields, and recalling that matrix elements of tensor operators 
vanish, the following operators can contribute to the B q — > r + r~ decay 



(O A IN y = 4GV (gD q (Qm Y b M ) (f N 7m r N ) , 
(CO* = 4GV {g^T {ft Y b L ) (q L 7m ti) , 
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(O^y = AG F (g^y (q M b N ) (f N r M ) , (M^N) 
(0{ L y = 4G F (g^ L r (q L Y h) (f L 7/l r L ) . (2.8) 

In these equations the subscripts 0, 1/2, and 1 denote the isospin of the field bilinears, and 
M, N = L,R label the fields' chirality. In terms of {gf IN ) q , the C% coefficients in Q2.5| ), which 
are directly constrained by the experimental data, read 

c q P = -(g$y - ( 9 ?/ L 2 y , cy = { g L x «y - ( g * L 2 y , 

C\ = {g L L y + ( 9 n g - (g^ R y - (g£ L y + (g^y + (g^ L y • (2.9) 

Some of the heavy states which generate the effective operators in (|2.8|) appear in non- 
trivial SU(2) multiplets. Assuming that the mass splittings between different members of 
the same multiplet are not too large, SU(2) rotations leave the overall coefficients in ( p.8| ) 
invariant to a good approximation. This allows us to obtain model independent relations be- 
tween contributions to different transitions of some of the operators in ( |2.8| ). The operators 
whose contributions to B — > r + r _ can be constrained in this way are listed in Table |. Op- 
erators corresponding to the effective couplings in the first column are related through SU(2) 
rotations to operators which induce the decays B — > X q v v, B — > r and B — > X q r u, as 
given in the second column. These decay modes provide the strongest constraints on various 
new physics contributions to B — > t + t~ . The bounds on BR(5 — > X q vv) (|1.3| ) and on 
BK(B — > t v) (|1.2|) imply that the contributions to B — > r + t~ proportional to the coeffi- 
cients in the first two lines of Table [I] are much below the present experimental sensitivity. 



Constrained operators 


Decay mode 


(9 A IL y, (tf L y 

(s^) 9 , (g^ LL ) q , (g^ L y 


B X q v v 

B^TU 

B — > X q t v 



TABLE I. Effective couplings of the operators that contribute to B q — > t + t , which are con- 
strained by the limits on the decays listed in the second column. 

6 



The coefficients {gi^) s and {gQ LL ) s in the third line are only weakly constrained by the data 
on BR(I? — > X c tu) ( [1.11 ). This is because the SM contribution to this decay is large, and 
its possible interference with the new physics cannot be neglected. Since the sign of the 



interference is not known, ( |1.1| ) does not yield definite limits on the model independent pa- 
rameters. Finally, the contributions proportional to {gQ IR ) q (M = L, R) and to (g^) 9 are 
presently unconstrained, as they are not related to any existing experimental limit. 

III. EXPERIMENTAL STATUS 

To date, no dedicated experimental search for B — ► r + r~ (X) decays has been carried 
out. In this section we discuss the possibilities of searching for these decays in current 
experiments. 

The CLEO Collaboration has established limits on decays into any pair of charged 



leptons |13[, including different final state flavors, except for the B& — > t + t~ decay mode. 
The reason is that in all but this case the final state contains a muon or electron with a 
well-defined energy that can be easily searched for. The CDF Collaboration has recently 



established strong limits on BR(l^,s — > /J> + fi~) || , These limits follow from the absence of 
muon pairs with invariant mass matching m# d or mg s . Because of these selection criteria, 
the CDF analysis does not constrain B r + r~ (X) followed by r — > \i. 

The stringent UA1 bound BR(I? — > X fi + fi~) < 5 x 10~ 5 |T5| has been used to constrain 
the product of branching ratios BR(i? — * X r + t~) x [BR(r — ► fiuu)] 2 , and thus to infer 
an indirect limit on BR(5 — > X r + t~) 0]. However, the UA1 Collaboration searched for 
muons pairs with large invariant mass, 3.9 GeV< m w < 4.4 GeV. Muons from r decays 



would not have passed this cut, so the limit inferred in [16| does not hold. 

We conclude that the existing data still allow for B — > t + t~ (X) branching ratios up 
to 0(10%). Therefore, in searching for B decays with a t + t~ pair in the final state, mea- 
surements at LEP can be competitive with other searches, and may even yield the best 
bounds until asymmetric B factories will start operating. Unlike at B factories running on 
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the T(4S), both Bd and B s meson decays can be studied at LEP. Since the LEP b hadron 
sample contains about 40% and 12% B s mesons ||17|| , the limits on B s decays are weaker 
than the limits on the corresponding Bd decays by about a factor of 3.3. 

The absence of B decays associated with very large missing energy, which yielded the 
limit on BR(5 — ► tv) (|1.2j ), also constrains BR(Z? — > r + r~). However, compared to the 
decay B — > r v, the B — > t + t~ mode yields a much softer missing energy spectrum, since 
in this case both neutrinos come from secondary decays. In addition, to reject background 
from semileptonic b and c decays, events with charged leptons in the final state are rejected. 
This weakens the limit by an additional factor of about 65%, corresponding to the hadronic 
t branching fraction. Still, for sufficiently large B —> t + t~ branching ratios some events 
would have been seen in the large E miss region studied by ALEPH || . To obtain the bound 
implied by the absence of such events, we need to evaluate the probability of B — > r + t~ 
decay events to pass the -E m i SS > 35 GeV cut, relative to that of B — > tv decays. We 
estimated this probability with a Monte Carlo simulation similar to that in 0, except that 
now we used the hadronic invariant mass spectrum in r decays as measured by CLEO|] We 
also made some simplifying approximations which are not always conservative, and could 
be avoided in a dedicated experimental analysis. For example, we neglected the effects of 
the correlation between the direction of the missing momenta from the r + and r~ decays. 
Nevertheless, we think that our results give a reasonable estimate of the limits that can be 
established by a dedicated experimental analysis of the LEP data. We obtain the following 
bounds: 

BR(5 d -> t+t-) < 1.5%, 

BR(5 S -> r+ t~) < 5.0% . (3.1) 

It is interesting to mention that the first of these limits is probably close to the bound that 
CLEO may be able to obtain using the fully reconstructed B decay sample ||18|| . As we shall 



^We thank Alan Weinstein for providing us with this spectrum. 
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discuss in the next section, in spite of being over four orders of magnitude above the SM 
predictions, the limits ( |3.1| ) yield the first constraints on some new physics parameters. 

Neutrinos from the B — > X r + t~ decay yield a missing energy spectrum which is too soft 
to produce any signal in the E m i SS > 35 GeV region. However, for large enough branching 
ratios, events from B->Ir + r would enhance the signal in the missing energy region used 
to measure the B — > X c t v decay. Taking into account that also for B — > X r + r~ selecting 
hadronic r decays weakens the limit by a factor of about 65%, we estimate that a bound 

BR{B ^ X t + t-) < 0(5%) (3.2) 

is within the reach of LEP sensitivity. 

Before concluding this section, we mention that the radiative decay B — > 7 v v is also 
associated with large missing energy. The corresponding branching ratios in the SM have 
been recently estimated to be of order 10 -9 for B& and 10~ 8 for B s [j!9 . We can obtain limits 
on these decays by assuming a missing energy spectrum for the B — > 7 v v decay, similar 
to that in B — > X v v in the limit of zero invariant mass for the final state hadron system. 
Taking into account that only neutral B mesons can decay into 7 v 9, while all b hadrons 
can decay via the b — > svv transition, we estimate the bounds 

BR(5 d ^7^)<lx KT 3 , 

BR(5 S -> 7 z/z/) < 3 x KT 3 . (3.3) 

Our estimates indicate that the rates for B — » r + t~ (X) decays can be constrained by 
the missing energy method only at the few percent level. By refining the details of the 
analysis (for example, by optimizing the -E m i SS cuts), dedicated experimental searches could 
probably establish better limits. However, regardless of possible such improvements, it is 
unlikely that the missing energy method could yield much stronger bounds. Therefore, it 
is appropriate to discuss whether similar (or better) limits can be obtained by different 
analyses. A back-of-an-envelope estimate shows that a limit competitive with our results 
could be established at LEP by looking for two charged leptons from r decays, coming from 



a secondary vertex in the hemisphere opposite to a tagged b. It seems to us that also the 
inclusive lepton spectrum in decays measured by CLEO |?0] cannot yield more severe 
constraints than those found above. Finally, the decay mode B — > 7 v v could be effectively 
searched for, by looking for the decay photon in the hemisphere opposite to a tagged b. Such 
a search may yield significantly better limits than our estimates ( |3.3p . To what extent some 
of these analyses could improve the constraints derived above can only be decided on the 
basis of more detailed experimental studies. 



IV. NEW PHYSICS 



In this section we study the constraints on new physics implied by the limits on 
B — > t + r~ (X) decays. By comparing ( |2.5|) with the limits on BR(Z? — > t + t~) given in 
7L{) we obtain the following constraint on the coefficients C P ', C q P and C q A : 



c q P + -c\ 



' cV 2 < 2.0 x 1(T 4 



+ 
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190MeV\ 2 /1.5ps 



Sb ) 



TB 



BR(B q 



r 1 r 



1.0 x 10" 2 



(4.1) 



3.0 xl0~ 4 ; for q = d, 
1.0 x 10~ 3 ; for q = s . 
In the next two subsections, we give two examples of models in which these bounds yield 
the first constraints on some parameters: models with light leptoquarks, and SUSY without 
R-parity. First we express in terms of the model parameters the effective couplings of the 
operators (gf 4N ) q in (|2.8| ), which arise from integrating out the new heavy states. Then, by 
inserting the expressions ( [2.9|) for the relevant C q coefficients into ( f4.1|) , we derive constraints 
on various couplings. 



A. Leptoquarks 

Leptoquarks (LQ) carry both baryon and lepton number, and hence couple directly 
leptons to quarks. A comprehensive analysis of the experimental constraints on the LQ 
couplings has been given in flIE|] , and is summarized in Table 15 of this reference. As 
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discussed above, the limit on BK(B — > X fi + ) [plj does not constrain BK(B — > X r + r ). 



Therefore, some of the limits on LQ Yukawa couplings in |Tb| involving the third generation 
do not apply. 

Several different types of LQ are possible, and most of them can induce the 
B — > t + t~ (X) decays. However, in many cases LQ also mediate the decays B — > X v v 
and B — > t u, which are tightly constrained. Therefore, we will concentrate only on those 
cases where transitions involving neutrinos are not induced. This can happen either be- 
cause of the particular electric charge of the LQ (for example \Q\ = 4/3), or when the LQ 
couplings to the left-handed lepton doublets vanish. We adopt here the notations of |16| . 



Scalar and vector LQ are denoted as S and V, while SU{2) singlets, doublets, and triplets 
are respectively labeled with a lower index 0, 1/2, and 1. The types of LQ relevant for 
B — > t + t~ (X) decays, and for which no strong constraints exist from other processes, are 

So, S 1/2 , Kf, Vf /2 . (4.2) 

The relevant scalar and vector terms in the interaction Lagrangian can be found in flIS 



Schematically, they read 

£ LQ = \-?£ l q J <pLQ, (4.3) 

where £i and qj denote respectively a lepton and quark fields, 4>lq represents one of the LQ 
in (|4.2|) , and i and j are generation indices. In deriving our constraints, we assume that 
all the \fj couplings are real, that only one type of LQ is present at a time (for other 
possibilities see |21]]), and we neglect the rotation from the interaction to the mass basis 
|22|| . Integrating out the LQ fields and Fiertz transforming, we obtain the expressions for 



the unconstrained coefficients of the relevant effective four-fermion operators given in (|2.8|) 



_ (^flK (Ar )33 (Ag%g (^r)33 

2 

So — 



(o R 

[9 ° ' 8G F m 2 ~ ' AGp m 2 

id") ' 



Lily, i^R /2 )3q (Afl /2 )33 (^R /2 )3q (^R /2 )33 

8G F ml, n ' AG F m\ 



°S 1/2 ™F i'o Vl/2 
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(9o 



iLL\q . 33 



AG F m 



Vo 



33 /. ,\ 

- — ^tt— 2 — ; \9i/2) - — or 2 — • 1 [ > 



For the different products of LQ couplings (A^) 32 (A^) 33 involving the second and third 
generation fermions, from (|2.9|) and ( |4.1|) we obtain the following limits: 

Aft A/? , A R Aft < 4.4 x 10 (j^^J , 



Ag Ag , A? A^° , A^ 1/2 A^ 1/2 < 2.2 x 1(T 2 (^g^) 



2 

j 

2 



A^° Ag , Ag A^° < 5.9 x 1(T 3 ( -^-,) , (4.5) 
l r > a l VlOOGeV/ ' v ; 

where the indices (32) and (33) respectively for the first and second coupling in the products 
are understood. The bounds on the analogous products (A^) 3 i (A^P) 33 involving a first 
generation (q = d) quark are about a factor of 1.8 stronger. From Table I we see that for 
q = d the products of the LQ couplings in ( f4.4p corresponding to (gQ LL ) 9 and (gf^) 9 are 
already tightly constrained by the limit on B — > r v decays, while for q = s they are weakly 
constrained by the upper bound on B — > X c t v. For all the other combinations, the bounds 
obtained from B — > r + r~ (X) are the strongest. 

B. SUSY without R-parity 

In SUSY models, it is usually assumed that R-parity is a good symmetry. However, this is 
not necessarily the case, and phenomenologically viable models have been constructed where 



R-parity is not imposed as an exact symmetry ||23|| . In the absence of R-parity, additional 
baryon and lepton number violating terms are allowed in the superpotential. Some of these 
terms can induce a large enhancement for certain rare B decay modes. Denoting by L l L , Q % L , 
£ l R , and d R respectively the chiral superfields of the z'th generation containing the left-handed 
lepton and quark doublets, the right-handed lepton, and the down-type quark singlet, the 
R-parity violating terms relevant for B decays read 
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Wft = X ijk U L V L I k R + \' l]k U L Q{ 4 . (4.6) 

The terms in ( [4.6|) give rise to two types of diagrams that can mediate B — > r + r~ (X). 
Exchanging a left- or right-handed squark gives rise to effective operators proportional to 
A' A'. Since squark exchange induces also the decay B — > X v u, these operators are already 
tightly constrained, and hence irrelevant for the present discussion. The exchange of left- 
handed sleptons generates operators proportional to A' A. These operators do not induce 
B — > X vv, but they can still contribute to B — > r v and to B — > X c r v. The corresponding 
effective couplings are: 

\l. n i A fc33 RL\q _ ^'k3q \ 

i% ' {9l ' 2) ~ 4G F m 

where k = 1, 2 due to the antisymmetry in the first two indices of the A couplings. Neglecting 
possible cancelations between the contributions from L\ and L 2 exchange, and assuming that 
one of the two couplings A'^ and \' kq3 is dominant, from ( |IO| ) and ( ^4.1| ) we obtain 



( „LR\q _ v fcg3 /v fc33 / RL\q _ yy k3q /y k33 i, 7 \ 



A' 23fc A fc33 , X' 32k A fc33 < 1.2 x 10- 2 ( M ^) ■ (4-8) 

The bounds on the analogous combinations involving the couplings of the d quark are about 
a factor of 1.8 stronger. From Table | we see that more stringent bounds already exist on the 
combination involving A' fel3 , while for A' fc23 additional weak constraints can be derived from 
B — > X c t v . However, for the products involving the X k3q couplings, these are the strongest 
limits. 



V. SUMMARY AND CONCLUSIONS 

To date no experiment has searched for B — > t + t~ (X) decays, thus no experimental 
bounds exist on the corresponding branching ratios. B decays into final states involving r 
leptons can be searched for by looking for the missing energy associated with the neutrinos 
from r decay. Such searches have been carried out at LEP to measure the B — > X c t v decay 
rate and to set bounds on the B —>■ r v and B — > X s v v decays 0,§ . In this paper 
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we pointed out that similar analyses can also set bounds on B — > r + t~ (X) decays. We 
estimated that limits at the few percent level [see eqs. ( |3.1| ) and ( |3.2| )1 are within the reach 
of the LEP sensitivity 

The SM predictions for the B — ► r + r~ (X) branching ratios are of order 10~ 6 or below. 
Therefore, the current experimental sensitivity only allows us to derive constraints on physics 
beyond the SM. To identify what kind of new physics could yield large enhancements, we 
performed a model independent analysis of these decays. While the decay B — > X r + t~ 
can depend on all the ten possible Dirac structures (12. If) , Bd, s — > t + t can be induced only 
by the three operators in ( |2.2|) . Thus, for studying new physics, these decay modes are 
complementary to one another. Much could be learned from the individual rates, and also 
from their ratio. 

Operators which can induce the B — > r + r~ (X) decays can be related to operators con- 
tributing to other processes through SU(2) rotations. In particular, when the dominant 
operators involve left-handed lepton doublets, B —>■ t + t~ (X) can be related to B — > X v v 
and to B — > tv. In these cases the limits on these processes 0,0] provide strong constraints. 
Therefore, only operators involving right-handed r's can induce B — > r + r~ (X) at the level 
of current experimental sensitivity. This restricts the types of new physics models that can 
be constrained through these decays. We gave two examples of such models: light lepto- 
quarks, and SUSY without R-parity. In both cases, the rather weak limits that we have 
estimated already yield the strongest constraints on some of the model parameters. 

In the future, B factories will provide significantly larger samples of B decays. Hopefully, 
the B — > r + r~ (X) decays will be observed, even if their rates are as small as predicted by 
the SM. Measurements of various B decay rates into final states involving r's (or iv's) a s well 
as of other observables (like the r polarization f24]| ) would provide very valuable information. 
Even if the experimental difficulty of these measurements will only allow establishing tighter 
limits on the decay rates, this will still yield strong constraints on possible physics beyond 
the SM, and might provide us with a better understanding of the third generation. 
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